Spin Relaxation in the Quantized Hall Regime 
in the Presence of Disorder 

S. Dickmann 

Institute for Solid State Physics of RAS, Chernogolovka, 142432 Moscow District, Russia 
and Weizmann Institute of Science, 76100 Rehovot, Israel 

We study the spin relaxation (SR) of a two-dimensional electron gas 
(2DEG) in the quantized Hall regime and discuss the role of spatial inhomo- 
geneity effects on the relaxation. The results are obtained for small filling 
factors {v <C 1) or when the filling factor is close to an integer. In either 
case SR times are essentially determined by a smooth random potential. For 
small u we predict a "magneto-confinement" resonance manifested in the en- 
hancement of the SR rate when the Zeeman energy is close to the spacing 
of confinement sublevels in the low-energy wing of the disorder-broadened 
Landau level. In the resonant region the -B-dependence of the SR time has a 
peculiar non-monotonic shape. If i/ ~ 2n+l, the SR is going non-exponentially. 
Under typical conditions the calculated SR times range from 10~® to 10^^ s. 

PACS numbers: 73.43.Cd, 72.25.Rb, 75.30.Ds 

1. For the relaxation of electron spins to occur, two conditions have to be fulfilled: 
the first one is the presence of an interaction mixing different spin states in the system 
studied; the second is the availability of a mechanism which makes the relaxation process 
irreversible. Both conditions can be realized in a rich variety of ways, and even in the case 
of two-dimensional (2D) electrons one finds a wide scatter of experimental data [1] and 
theoretical results [2-10] devoted to SR problems. Besides, in the magnetic field the SR 
process is actually the relaxation of the Zeeman energy \gfiBBASz\ {B \\ z, AS^ = — Sq 
is the spin deviation of the Sz component from the equilibrium value So). 

Theoretically, the relaxation problem of a flipped spin in a semiconductor heterostructure 
in high perpendicular magnetic field seems to be first formulated by D. Frenkel [4], who 
considered the relativistic part of the phonon field acting directly on the electron spin as the 
interaction mixing the electron spin states. One could estimate that in the case of another 



mechanism, exactly due to the spin-orbit (SO) couphng reformulated for the 2D case [2,3], 
the relevant spin-flip transition matrix element is at least by an order of magnitude greater. 
However, the work of Ref. [4] caused some misunderstandings even almost a decade after its 
publication. (See the comment [9] and references therein). 

Any properties of a 2DEG in the quantum Hall regime crucially depend on the filling 
factor V = N/Nff), where N and A^^ = L^/27r/^ are the numbers of electrons and magnetic 
flux quanta, respectively (L x L is the 2DEG area, Ib is the magnetic length). In this Letter 
we consider the SR problem in two formulations. 

First, we solve the problem of one- electron spin-flip in the presence of a random potential, 
and in so doing we study the SR in a lateral quantum dot in high magnetic field, where the 
effects of interplay between the locahzation and the Zeeman couphng are essential. We 
emphasize that we deal with a weak confinement but in the presence of strong B. Therefore 
well smaller energies are relevant (~ 1 K) than, e.g., in the case of usual interplay of Fock- 
Darwin states [11]. (The latter is based on competition of confinement and cyclotron energies 
which are >;10K). The most striking manifestation of the "magneto-confinement" effect 
occurs when the Zeeman energy is close to the spacing of lowest-energy levels in a quantum 
dot or in a potential minimum. This results in strong enhancement of the SR rate. The 
studied one-electron problem is actual for small filling factors, u <^ 1, but in strong magnetic 
fields it may also be extended in terms of the Hartree-Fock approach to a vicinity of an even 
filling factor, when u ~ 2n. 

Second, we report on the results for the SR in a strongly correlated 2DEG when it presents 
a quantum Hall ferromagnet (QHF), i.e. when the filling factor u is close to an odd integer 
[v ~ 2n In this case we study a new mechanism of relaxation of the total 2DEG 

spin, where the cause of irreversibility is neither electron-phonon interaction [6] nor inter- 
spin-wave scattering [7] but a disorder. We will see that in real experimental regimes this 
disorder relaxation channel has significantly to prevail over the phonon one. 

In both SR problems the temperature is assumed to be equal to zero which actually 
means that it is lower than the Zeeman energy. 

2. We consider a smooth random potential (SRP) as the disorder. The total single- 
electron Hamiltonian is thereby as follows: !K — ^^q^/2m* — ez&z/'^ + ^(r) -|- Hso + t4-ph, 
where ez = \g\iiBB, q = — iV -|- eA/c and r = (x, y) are 2D vectors, ^(r) is the SRP field. 



the Hso and C/g-ph terms respond to the SO and electron- phonon interactions (see below) . If 
the SRP is assumed to be Gaussian, then it is defined by the correlator K{v) — (w(r)ti(O)). 
We choose also {u{v)) = which means that the SRP energy is measured from the center 
of the Landau level. In terms of the correlation length A and Landau level width A, the 
correlator is 

i^(r) = A2exp(-rVA'). (1) 

In the realistic case A ~ lOK, A ~ 30 — 50 nm, therefore A>^^/m*A^. We study the case 
A <^ hide i^c is the cyclotron frequency), and Ib- In the SRP field the electron drifts 
quasiclassically along an equipotential line. However before the spin-flip it relaxes to a SRP 
minimum. Estimates for this relaxation time (due to phonon emission without any spin-flip) 
yield the values not exceeding 1 ns. 

For simplicity, we model the SRP in the vicinity of a minimum by a parabolic conflnement 
potential u — m*a;^r^/2, and to describe the electron states we use the symmetric gauge basis 
(A = rxB/2): 



(L^ is a Laguerre polynomial; note also that only the states with n+m > are considered 
in the following.) For the length a it should be substituted a — (/i/2m*Q)-'^/^, where Q — 



^/u'^+uj^/A. The system thus becomes equivalent to a lateral quantum dot [8,10], and 
we deal with the Fock-Darwin states [11] with energies En^m — h{2n + m + l)Q — hWc'm/2. 
The appropriate quantity is also the level spacing 5 — h{fl — u!c/2) ^ huj'^/ujc, concerning 
the m > states which belong to the same number n. We calculate the total rate of the 
transition of an electron initially occupying the upper spin sublevel to any final state of the 
lower spin sublevel. At flrst sight we should consider the spin-flipped state |0,0, J,) as the 
initial one. However, a correction of the states due to the SO coupling has to be taken 
into account. We use the SO Hamiltonian specifled for the (001) GaAs plane: Hso — 
q; (q X a\ + (5 {qy^y — qx^x)- This expression is a combination of the Rashba term [2] (with 
the coefficient a) and the crystalline anisotropy term [3,5-10] {cF^^y^z are the Pauli matrices). 
Assuming that a and (5 are small {a, (3 <^ fiuJclB) we flnd after perturbative treatment the 
spin-orbitally corrected states before and after the spin flip: 

N) = C'i,i|0,0,i) + C2,i|0,l,T) + ^^^y|^|l,-l,T) and (3) 



|/n.) = C'i,^|0,mJ) + C'2,^|0,m-14)-— ^|l,m-14) + -^^ (4) 

where < m < e^/^ (the Zeeman energy ez has been neglected as compared to hio^. The 
coefficients Ci^m are defined as follows: let T = ab^j ah£L{ez—S) and Pm — 2\/l-\-T'^m, then 
Ci^m = \/l/2 + l/Pm, and C2,m = i sign {S — ez)\/l/2 — l/ P^- Here the resonance mixing 
of the "spin-up" and "spin-down" states (if ez — S) has been properly taken into account. 
Note that the behaviour of the states (3)-(4) in the vicinity of the resonance (namely in the 



interval AB/B ^ ayml/h^ujc ^ 0.1) is governed only by the Rashba SO mechanism. 

In the resonance region Ci^i ~ |C'2,i|, and \i) is thereby a well hybridized spin state. The 
~^ l/o) transition then is not due to the SO couphng, and this should lead to the SR 
enhancement. The final state |/o) is the almost "pure" spin-state. In fact we may always 
set Ci^rn — 1 and C2,m = in the expression for \fm)- Indeed, though the spin hybridization 
of the state is significant, it plays negligible role in the SR process because of vanishing 
of the relevant phonon momentum {ez— S)/cs. Then we find the matrix element (/m|^e-phK) 
and, using the Fermi golden rule, obtain the SR rate within a certain SRP minimum, I/t^^ — 
fj:m,Mm\U^ph\i)f5ihcsk - ez). Here C/e.ph(r) = (Vy)V2 [/^(q, A;,)e^q^ where the 
index s labels phonon polarization, V is the sample volume, and Ug is the renormalizcd (in 
the 2D layer) vertex which includes the deformation and piezoelectric fields created by the 
phonon [6,12]. The summation over s involves averaging over directions of the polarization 
unit vector for both components of the clcctron-phonon interaction, and this may be reduced 
to IE. Usf - nhcsk/plTAik), where r^"^ = r^' + ^T-^pl{q^kl + qUD/k'' (see Ref. [6]). The 
nominal times for the deformation and piezoelectric interactions in GaAs are td ~ 0.8 ps, 
and Tp ~ 35 ps [6,12]. The nominal momentum is po = 2.52 ■ 10^/cm [12]. (We also refer 
to Refs. [6,12] for details concerning the meaning of these quantities and their expressions 
in terms of the GaAs material parameters.) Finally, in the general expression for I/t^j we 
perform the summation and after a routine treatment arrive at the result 

u2 



- = ^ / ^ E ^ [^UO+^m (5n. + 5^-35^/8)6-^"^ (^-Pl , (5) 
Sapor, 7o VWo<.^^/5^' V 2 ; 

ez-mS bljp f 2m /7, \ i , "^m fl ■ ^ 

where bm^a— , = Am= ^-^/^bra C2,i y 77 sign (d-e^), 

hcs [apoyrd \bmV^ ) '^^^ V 2 

and ^m — T-pr-\ - ( ^~Ci,iT~' — ) • '^^^ estimate for a and (3 depends on the effective layer 



hQa\2\ ^'^S+ez, 
width [2,3]. The rate 1/r^ as a function of 5 at a; = 5K is shown in the inset of Fig.l for 



realistic parameters indicated in the capture. 

At a; = (i.e. in the "clean" limit), the summation in the formula (5) is carried out over 
all numbers m = 0, 1, ...oo, and the expression (5) is reduced to 

- - V /' -^e-'-^/^S + 5C - 35CV8) , (6) 

To Jo V 1 - ? 

where V = ^-e— oc B, b— — — oc VB, S = — oc B . The bold curves m 

ApohrujlclTp ncs td [ncspo ) ^ 

the inset and in the main picture of Fig. 1 show the corresponding S-dependencies. 

Note that if ez = 0, then at any tq the projection V^tq) = (/m|(5(r— ro)|i) vanishes when 
calculated at a finite cu in the leading order in the SO constants. [It does not occur for the 
"clean" states, i.e. if in Eqs. (3) and (4) we pass to the a; ^ hmit before equating €z 
to zero.] Such a vanishing is a manifestation of the general feature [8,10]: at zero Zeeman 
energy the effects of the SO coupling in the leading order and of the orbital magnetic field are 
similar. In particular, in a quantum dot (where 5 > e^!), the first order SO approximation 
in the ez ^ limit results only in a small rotation of eigen states in the spin space. With 
decreasing B we get V oc B^^"^ (if "C 5 <^ <^c) and obtain a sharper fall of the relaxation 
rate as compared to the "clean" case (6). 

In the presence of the SRP we have to carry out averaging I/tcs = dLjF{Lj)/T^, 
where the distribution function F{uj) is the probability for the confinement frequency to 
take a certain value cu. One may prove that in the case of a Gaussian potential u{r) it should 
be chosen in the form F{u)) = -^^^^ cxp {—uj'^ / AlOq) . [The value co^ is proportional to the 
curvature V^m, and a routine analysis yields (jJq = 6^/''(A/m*)^/^/A.] Calculating 1/Tes with 
this function, we obtain the final result (see Fig.l). As it has to be, in comparison with r^^ 
the resonant behaviour of l/reff is smoothed, but at actual values of cuq it results in non- 
monotonic 5-dependence of Teff. Beyond the resonance region the behaviour is as follows: 
(i) at small magnetic fields (when ez <^ HcUq/cUc <^ cUc) only one final state |/o) participates 
in the SR, and we find that Teff oc B'^; (ii) at high fields (when ez ^ Sq) there is a large 
but finite number of possible states \fm) into which the confined spin- flipped electrons could 
relax, and the SR time is always longer than tq but approaches this with increasing magnetic 
field. Note that exactly in this high-field regime the one-electron model becomes relevant for 
fillings v ~ 2n. Then the total 2DEG spin is determined only by a small amount \p—2n\N(f, 
of effectively "free" electrons/holes belonging to the (n-|-l)-st/n-th Landau level. 



3. So, the problem has been solved when the total 2DEG spin is well smaller than 
N^. Now we study the opposite case: in the ground QHF state the spin numbers attain 
maximum, S = Sz = This case is also remarkable, since to the first order in the ratio 

Tc — {e^ /elB)/fiu!c the low- lying excitation are again exactly known: these are 2D spin waves 
or spin excitons (SEs). The most adequate description of the SE states is realized by spin- 
exciton creation Ql^ = e~*^''^'Ba|^2£'^tp-^ annihilation Qq = (Qq)^ operators 

[13]. In this definition a^p stands for the Fermi annihilation operator corresponding to the 
\n,p, a) — L~^/'^ ^pn{x-\-piy) state in the Landau gauge is the nth harmonic oscilatory 
function). So the one-exciton state is QjjiO), where |0) stands for the ground state. At small 
2D momentum [qIb -C 1) the one-exciton state has the energy 8q — ez+^qlBY /'^^n (now 
we need only this small momenta approximation, see also general expressions for the 2D 
magneto-excitons at integer fiUing factors in Refs. [14]). M„ is the SE mass aX v = 2n + 1, 
namely in the r^^oo limit: 1/Mq = (e7£^B)vW8> l/^i ^ T/^Mq,... Note that the sum 

— Ylii^^- lowering the spin number by 1 {i labels the electrons), when considered 
to be projected onto the nth Landau level, is simply proportional to the "zero" exciton 
creation operator Q\. When the SO coupling is ignored, any state (S'_)^|0) is the eigenstate 
independently of the Tc magnitude and of the presence of a disorder. 

The SO interaction Hso and the SRP field u{r) may be accounted perturbatively as 
usual. Meanwhile, now the unperturbed part of the Hamiltonian involves the Coulomb 
interaction. For present purposes the approximation in the context of the projection onto 
a single Landau level is quite sufficient (c.f. Ref. [6]). We note also that the QHF state 
|0) is resistant to the SRP disorder. Such a stability is determined by the exchange energy 
(~ e^/ kIb) which is much larger than the amplitude A. 

There are two fundamentally different alternatives to provide the initial perturbation of 
spins. The first one is the perturbation of the spin system as a whole when the S number 
is not changed: AS" = 0, but A5'^ ^ 0. This is a Goldstone mode which presents a quantum 
precession of the vector S around the B direction. In terms of the SE representation | A5'^| = 
N^/2—Sz is the number of "zero" SEs excited in a 2DEG. Let lA^'^l =N] the corresponding 
state is {SJ) |0) oc (Qo) |0), and it has the energy Nez- [Note that "zero" SEs do not 
interact among themselves; besides, the stabihty of the (5'_)^|0) state with respect to the 
disorder is identical to that for the |0) state because of ^(r) and S_ commuting.] The second 



case of the perturbation is the AS — AS^ type of the deviation. This does not change the 
symmetry and involves the excitation of "nonzero" SEs, where each SE changes the spin 
numbers by 1: S Sz~^ Sz—l- In contrast to '"zero" SEs, the ''nonzero" ones interact. 
This interaction [7] or/and the direct exciton-phonon couphng [6] govern the "nonzero" SE 
annihilation which should go faster than for "nonzero" SE annihilation process. 

The SRP inhomogeneity does not essentially affect the SE energy and the "nonzero" 
annihilation. Indeed, the exciton is neutral, and the interaction with the SRP incorporates 
the energy C/x-srp ~ qI'b^I ^ (^^^ "nonzero" SE possesses the dipole momentum eZ^[qxi], 
see Ref. [14]). The latter may be negative but is always smaller than the "nonzero" SE 
energy £q. If the SE would annihilate, the energy conservation condition £q + t/x-sRp — 
is not satisfied except of negligibly rare points, where the gradient Vu is accidentally large. 
Therefore, the SRP leads only to small corrections (of the order of C4-sRp/^^g) to the other 
mechanisms of the "nonzero" SE relaxation [6,7]. 

A distinctly different process contributes to the SR in the case of the first type of the spin 
perturbation. This is an effective interaction of "zero" SEs among themselves arising due 
to the SO coupling [6]. Such an interaction does not preserve the total number of excitons 
N: at elementary event two "zero" excitons merge into one "nonzero" (the spin momenta 
change following the rule + S^S — \). In other words, the SO coupling and the 

SRP field mix the initial \i) = (Qo) |0) the final \ fq) = \Qq) |0) states. These 
"many-exciton" states have to be normahzed (see the normalization factors in Refs. [6,13]). 
So, the SRP plays the same role as the phonon field studied previously [6]. Now the energy 
conservation condition takes the form: 2ez ~ Sq, where the interaction of "nonzero" SE with 
the SRP is ignored as compared to other members of this equation. If solving this for g, we 
obtain q = qq = \/2Mnez/lB- 

The detailed calculation of the SR rate is truly similar to that performed in Ref. [6]. 
Indeed, the Fourier expansion u{r) — J^q^io)^^^'^ looks like the phonon field created by 
"frozen" (of zero frequency) phonons. Then this SRP field and the SO Hamiltonian are 
treated perturbatively. In so doing, it is convenient to present them in terms of the excitonic 
operators. We calculate the relevant matrix element between the \i) and \fg) normalized 
states and obtain |Mi^/J =N^{a'^ + (3^)\qu{q)\^/{huJcyN^. Finally, again with the use of 
the Fermi golden rule, we find after summation over all states \fq) that the SR rate takes 



the same form as in the case of the phonon mechanism [6], dSz/dt — {^Szf' /ts^vN^^ but 
incorporates a different time constant: 

.1 l&7:\a^ + l3^)MlezK{qo) 



Here K stands for the Fourier component of the correlator [the equivalence K[q) — 
LF' \u{q)f' /At:'^ has been employed]. The SR follows the low 

where = TpJ^^+Tg~p, because the relaxations of both types proceed in parallel. A natural 
question is: what is the ratio of the times Tph and r^nv ? K T^hCsQa (in fact, this means that 
T < 1 K) and B <15, then the SR time Tph depends weakly on T and B. In particular, at 
1^ = 1 we find that Tph — 10 /is [6,7]. The ratio of interest is determined only by the Fourier 
component Kiqo): Tph/rgRP = OMiTTpPoA'^A^-^^io/^/h'^Cs [for K we have substituted the 
value calculated with help of Eq. (1)]. So, for the actual parameters Tph/rgRp ~ 100 — 1000, 
i.e. exactly the "disorder" time (rsRP~10~^ — 10~''s) governs the breakdown of this Goldstone 
mode. In conclusion, we remind that the relaxation is going non- exponentially and the actual 
time is increased by a factor of Nfj^/ ASz{0). 

4. I thank Y.B. Levinson for helpful discussion and acknowledge support by the MIN- 
ERVA Foundation and by the Russian Foundation for Basic Research. I thank also the 
Weizmann Institute of Science (Rehovot) for hospitahty. 

[1] M. Dobers, K.v. Klitzing, J. Sneider et al, Phys. Rev. Lett. 61, 1650 (1988); S. Bar-Ad, and 
I. Bar-Joseph, Phys. Rev. Lett. 68, 349 (1992); V. Srinivas, Y.J. Chen, and C.E.C. Wood, 
Phys. Rev. B 47, 10907 (1993); V.E. Zhitomirskii, V.E. Kirpichev, A.I. Filin et ai, JETP 
Lett. 58, 439 (1993); B.W. Alphenaar, H.O. Miiller, and K. Tsukagoshi, Phys. Rev. Lett. 81, 
5628 (1998). 

[2] Yu.A. Bychkov and E.I. Rashba, JETP Lett. 39, 78 (1984). 

[3] M.I. D'yakonov and V.Yu. Kachorovskii, Sov. Phys. Semicond. 20, 110 (1986)]. 

[4] D.M. Prenkel, Phys. Rev. B 43, 14228 (1991). 

[5] G. Bastard, Phys. Rev. B 46, 4253 (1992). 

[6] S.M. Dikman and S.V. lordanskii, JETP 83, 128 (1996). 



[7] S.M. Dikman and S.V. lordanskii, JETP Lett. 70, (1999). 
[8] A.V. Khaetskii and Y.V. Nazarov, Phys. Rev. B 61, 12639 (2000). 
[9] A.V. Khaetskii, Phys. Rev. Lett. 87, 049701 (2001). 
[10] LL. Aleiner and V.L Fal'ko, Phys. Rev. Lett. 87, 256801 (2001). 

[11] V. Fock, Z. Phys. 47, 446 (1928); C.G. Darwin, Proc. Cambridge Philos. Soc. 27, 86 (1930). 

[12] V.F. Gantmakher and Y.B. Levinson, Carrier Scattaring in Metals and Semiconductors (North- 
Holland, Amsterdam, 1987). 

[13] The excitonic operators appeared originaUy in pubhcations of A.B. Dzyubenko and Yu.E. 
Lozovik, Sov. Phys. SoUd State 25, 874 (1983) [ibid. 26, 938 (1984)]; in the paper of J. Phys. 
A 24, 415 (1991) these authors considered the many-exciton states and calculated the norms. 

[14] LV. Lerner and Yu.E. Lozovik, Sov. Phys. JETP 53, 763 (1981)]; Yu.A. Bychkov, S.V. lor- 
danskii, and G.M. EHashberg, JETP Lett. 33, 143 (1981); C. Kalhn and B.I. Halperin, Phys. 
Rev. B 30, 5655 (1984). 




Magnetic field, B (T) 

FIG. 1. Calculations are carried out for a = /?/3 = 10~^K-cm and Cg = 3.37 • lO^cm/s (other 
material parameters are given in the text). In the inset the position of the SR peak corresponds 
to the condition ez = 6. The evolution of i?-dependences of the spin relaxation time TeS with the 
parameter uq is shown in the main picture. 



